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^j. The statistical properties of coherent radiation scattered from phase-ordering materials are stud- 

' ied in detail using large-scale computer simulations and analytic arguments. Specifically, we consider 

a two-dimensional model with a nonconserved, scalar order parameter (Model A), quenched through 
an order-disorder transition into the two-phase regime. For such systems it is well established that 
the standard scaling hypothesis applies, consequently the average scattering intensity at wavevector 
k and time r is proportional to a scaling function which depends only on a rescaled time, t ~ |k| 2 r. 
We find that the simulated intensities are exponentially distributed, and the time-dependent aver- 
age is well approximated using a scaling function due to Ohta, Jasnow, and Kawasaki. Considering 
fluctuations around the average behavior, we find that the covariance of the scattering intensity for 
a single wavevector at two different times is proportional to a scaling function with natural variables 
St = \ti — ta| and i = (ti + ^2)/2. In the asymptotic large-i limit this scaling function depends only 
on z = St/i 1 / 2 . For small values of z, the scaling function is quadratic, corresponding to highly per- 
£jj ■ sistent behavior of the intensity fluctuations. We empirically establish that the intensity covariance 

(for k 0) equals the square of the spatial Fourier transform of the two-time, two-point correlation 
function of the order parameter. This connection allows sensitive testing, either experimental or 
numerical, of existing theories for two-time correlations in systems undergoing order-disorder phase 
I . transitions. Comparison between theoretical scaling functions and our numerical results requires no 

' O ■ adjustable parameters. 
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A scattering experiment, using neutrons or X-rays for example, is one of the most direct measures of the structure 
of materials. Naively, this comes about because in the Born approximation, which usually applies for X-rays and 
neutrons, the intensity in scattering measurements is proportional to the Fourier transform of a density-density 
[~^. ' correlation function. It is the wavelike properties of the scattering probe which produces the Fourier transform. For a 
OS . deeper understanding of the relationship between scattering intensity and structure one must realize that this direct 
correspondence applies precisely only for coherent waves. Indeed, for conventional sources, a given point in the incident 
wave is only coherent within a small volume of neighboring points. This coherence volume has transverse dimensions 
determined by how parallel the wavefronts are and a longitudinal length determined by how monochromatic the wave 
is. In a standard scattering experiment the different coherence regions of the incident beam scatter independently. 
The intensity measured thus depends on an incoherent average over different regions of the scattering volume. By 
restricting the scattering volume of the sample to less than the coherence volume of the beam, one can eliminate this 
Q , incoherent average, and thus learn more about the material's structure. Of course, the experimental difficulty which 
arises is to obtain sufficient diffracted intensity to measure a signal. Recently, it has been demonstrated that coherent 
diffraction experiments can be performed with X-rays using high-brilliance synchrotron sources Q. For coherent 
diffraction, the scattering from an inhomogeneous material displays a characteristic speckled scattering pattern. For 
instance, the random distribution of phase-ordering domains shown in Fig. 1 and discussed below, results in the 
speckle pattern shown in Fig. 2. As the domains change shape, the speckle pattern changes, and this time dependence 
of the speckle offers a unique method for studying the evolution of inhomogeneous materials. 

Motivated by these advances in experimental techniques, we have undertaken a theoretical study of intensity 
fluctuations caused by scattering from a nonequilibrium system undergoing phase ordering by domain growth. When 
a disordered homogeneous material is rapidly brought to a new set of conditions, corresponding to the coexistence of 
two equilibrium phases, a spatial pattern of domains of the two phases develops ||. This change of conditions is often 
accomplished by a rapid quench from a high temperature to a low one below a miscibility gap. The quench from a 
homogeneous state with local fluctuations creates a microstructure of interconnecting, interlocking domains through 
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the kinetics of a first-order phase transition. As time goes on the domains grow, so as to minimize the area of the 
domain walls that separate the phases. 

When the average domain size R(t) at time r is large compared to all other relevant lengths, except for the extent 
L of the system itself, the system looks invariant if all lengths are measured in units of that domain size. In this case 
the structure of these many domains is said to "scale" with R(t). Experimentally, the growing domain structure is 
often studied by means of the scattering intensity whose width is proportional to the inverse of R(t). When the 
average scattering intensity is scaled in units of this time-dependent length, one obtains the scattering function for late 
times in the form of a time- independent scaling function. The time-dependence then enters only through i?(r), which 
can typically be described in terms of an exponent n, such that R ~ r™. This scaling hypothesis has been found to 
apply to a large range of systems, and to be unaffected by many of the microscopic details of specific materials. That 
is, the scaling function and the growth exponent n are two features which are common to a large number of systems, 
collectively called a universality class. Universality classes for phase ordering by domain growth are delineated chiefly 
by the presence or absence of conservation laws. Indeed, many aspects of this nonequilibrium process can be described 
by relatively simple theoretical models. For example, for systems described by a nonconserved scalar order parameter, 
often called Model A ||, the growth exponent is found to be n = 1/2. Systems included in this class are the Ising 
model with spin-flip dynamics, binary alloys undergoing an order-disorder transition, and some magnetic materials 
with uniaxial anisotropy. Model B |2j refers to systems in which the scalar order parameter is conserved, and the 
only growth mechanism is diffusion. Systems in this universality class have n = 1/3 and include the conserved Ising 
model, as well as binary alloys undergoing phase separation. 

In the present paper we investigate the time-dependent fluctuations around this scaling behavior, and we demon- 
strate how to study these fluctuations experimentally through analysis of the time-dependent scattering. The first 
theoretical study of such behavior is given in Ref. p|. The scattering intensity is related to the Fourier transform 
of the order parameter ip(r, r), the scalar field describing the inhomogeneity of a specific sample of the scattering 
material, by 

J(k,T) = |^(k,r)| 2 , (1.1) 

where we ignore the proportionality constant for convenience. The average of /(k, r) over an ensemble of initial 
conditions is the structure factor, 

S(k,r) = (I(k,r)) . (1.2) 

Here, the ensemble average expresses the distinction between coherent scattering, given by /, and incoherent scattering, 
given by S. Fluctuations around the average scattering are the main topic of interest in this paper. The structure 
factor can also be expressed as the Fourier transform of the correlation function of ip(r, r), 

C(r,T) = <V(0,r)V(r,r)). (1.3) 

Because of these relationships, the scattering intensity and structure factor have been important tools for studying 
the dynamics of materials far from equilibrium. If a system is isotropic, then its average scattering properties depend 
only on the magnitude of the scattering wavevector, |k| = k. 

As mentioned above, scaling by the domain size R ~ r™ implies that correlations are time independent when 
measured in units of R. Specifically, for an isotropic system, C(r, r) = C(r/R(r)), where C(r) is one form of the 
scaling function. Fourier transformation gives the average scattering intensity in terms of another form of the scaling 
function, which depends only on a scaled time t oc rfc 1 /™, 

k d S{k,r) =F x {t) , (1.4) 

where d is the spatial dimension of the scattering material. The specific form of Fi(t) depends on the dynamic 
universality class. 

As discussed above, when the scattering involves coherent radiation, /(k, r) is directly measured without self- 
averaging. The scattering from an inhomogeneous material displays a characteristic speckled scattering pattern much 
like the one in Fig. 2, which shows the squared norm of the Fourier transform of the configuration in Fig. 1. The 
intensity of each speckle in /(k, r) is due to correlations in the inhomogeneous sample, and I — (I) = I — S gives the 
fluctuations around the average scattering intensity. 

For materials in equilibrium, deviations of / around the structure factor are induced by thermal fluctuations in tp. 
In fact, fluctuations in the speckle patterns from scattered laser light have been used as the standard basis of photon 
correlation experiments at wavelengths ranging from the ultraviolet to the infrared |5|||. With the advent of high- 
brilliance synchrotron photon sources, experiments involving coherent X-rays have become possible. For example, 
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Brownian diffusion rates in gold colloids have been determined from the time required for a change in an X-ray 
speckle pattern J?]]^] . Speckle from coherent X-rays has also been used to study equilibrium fluctuations in Fe3 Al near 
an order-disorder transition and in micellar block-copolymer systems |[o]]. X-rays can probe materials on much 
smaller length scales than is currently possible with lasers, and their greater penetration allows the study of optically 
opaque materials. 

In the present paper, we present a theoretical study of fluctuations in the scattering intensity from a nonequilibrium 
system undergoing phase ordering by domain growth. For such systems the intensity fluctuates around a time- 
dependent structure factor |4| Jll] , p^| . The time evolution of the intensity at a specific wavevector for one particular 
quench can be normalized by the average behavior. A typical example of such a normalized time series, obtained 
from a simulation at zero temperature, is presented in Fig. 3. It is the correlations of such time series, averaged over 
many individual speckles and quenches, that can be used to study the pattern formation process in phase-ordering 
materials. 

In Fig. 3 we also show a "Brownian" function which was constructed to have the same single-time probability density 
and an exponential two-time covariance with the same characteristic time as the normalized nonequilibrium scattering 
intensity. Qualitative differences are immediately evident in the two time series. The Brownian function fluctuates 
quickly, with large amplitude variations, while the intensity fluctuations produced by the phase-ordering system 
vary slowly, with markedly less variation in amplitude on short time scales |l2| . This property of the nonequilibrium 
intensity fluctuations is called persistence |l3| , and it indicates a qualitative difference between the two-time correlation 
functions for the two processes Jl4| . 

For the remainder of this paper, we specialize to Model A as a simple model for systems undergoing phase ordering 
following a quench through a second-order order-disorder phase transition. Some preliminary numerical results from 
simulations less extensive than the ones used here were presented in Ref. |l5[ . 

The outline of the rest of this paper is as follows. Section || describes the details of our numerical approach, which 
uses a standard time-dependent Ginzburg-Landau equation with a nonconserved order parameter. In Sec. [II the 
time-time covariance of individual speckle intensities is discussed, using analysis readily adaptable to experimental 
data. The equality between this covariance and the square of the two-time structure factor of the nonequilibrium 
material is argued in Sec. IV, and comparisons to specific two-time theories follow in Sec. There an analytic 
expression for the universal scaling form for the intensity covariance at late times is obtained. Finally, our results are 



summarized in Sec. VI 



II. METHOD 



The dynamics of speckle in k-space were simulated by generating successive scattering patterns from a real-space 
simulation of the dynamics of phase ordering following a quench through an order-disorder phase transition. The 
configuration of the real-space system is described by a nonconserved scalar order-parameter field ip(r,T), and its 
dynamics are governed by the following time-dependent Ginzburg-Landau equation, 



, SF[iI>(t,t)] 



+ C(r,r) 



(2.1) 



The first term on the right-hand side of this Langevin equation corresponds to deterministic relaxation, with rate 
constant T, towards a minimum value of the free-energy functional T[ip{v,T)\. Thermal noise, which we neglect, is 
modeled by the random variable £, whose intensity is proportional to T and temperature by virtue of a fluctuation- 
dissipation theorem We neglect £ because the most important sources of noise here are the initial conditions, 
which give the random domain morphology. The main effect of thermal noise at low temperatures far below any 
critical point, is only to thermally roughen the domain walls. 

We employ the standard Ginzburg-Landau- Wilson free energy 0j , 



dr 



y(r,r) + ^ 4 (r,r) + ^|V^r,r)| 2 



(2.2) 



For a > 0, the local part of the integrand in Eq. ( |2.2|) represents a bistable potential, and the parameters of the model 
define an equilibrium field magnitude, |^>o| — \Z a /u, a thermal correlation length, £o — \/2c/a, and a characteristic 
evolution time, r = (ar) _1 . The model can be rescaled by choosing a new field tp = tp/\ip \, a new length f = -y/^r/^o, 
and a new time f = t/tq. Dropping the tilde from the rescaled quantities for convenience, we get the rescaled 
dynamical equation, 
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^l^ = (l + V 2 )^(r,r)-^(r,r), (2.3) 

which has no adjustable parameters. Since thermal fluctuations are explicitly omitted in the simulation, the only 
randomness comes from the high-temperature state the system is in before the quench. This was implemented by an 
initial condition such that ip(r,0) consists of independent random numbers uniformly distributed between ±0.1. 

The simulations were conducted on square lattices with periodic bou ndary conditions, lattice constant Ar = 1 , 
and a system size of L x = L y = L = 1024. The Laplacian in Eq. ( |2.3| ) was implemented using the eight-neighbor 
discretization |l6|]l7| 

= 2^)2 ( E + lYl ^ nnn - > ( 2 - 4 ) 

where V'nn are the four nearest neighbors of ip (along the lattice directions), and V'nnn are its four next-nearest 
neighbors (diagonally). A simple Euler integration scheme with Ar = 0.05 was used to collect data up to a maximum 
rescaled time of r = 2000 for 100 separate sets of initial conditions. 

Usually, the order parameter takes the values ±|i/ , o| everywhere except at the domain walls, which are negligibly thin 
compared to the domains themselves. After our rescaling, domain walls have a soft nonzero width of approximately 
y/2. To minimize the effect of this nonzero, though small, width, we use a nonlinear mapping of the rescaled order 
parameter to ±1 before taking the Fourier transform. The transformed field, tjj, is defined by 

^(k, r) = -L= £ sign (V(r, r)) e lkr , (2.5) 



L d 

where the fact that the lattice spacing is unity in all directions has been used. The Brillouin zone in two dimensions is 
defined by the discrete set of wavevectors, k x , k y — 2nj/L with j e {0, ±1, ±2, . . . , ±{L/2 — 1), L/2}. The scattering 

intensity J(k, r) = |j/>(k, r)| 2 and its average over initial conditions is the time dependent structure factor, S(k,t), 
as described in Sec. Q. To be consistent with the numerical integration, the magnitude of the scattering wavevector, 
fc(k), is defined using the operator relation 



fc 2 (k)^(k,r) 



-l=^ e ^V 2 ^r,r). (2.6) 



Substituting the discrete version of the Laplacian from Eq. (2.4), one obtains for d = 2 



fc 2 (k) = 3 — cosk x — cos k y — — cos (k x + k y ) — — cos (k x — k y ) . (2.7) 

This method has been used to calculate the magnitude of the wavevector for scaling structure factors and data binning. 
Because of lattice effects we consider only those wavevectors with < fc(k) < 0.75. 

In order to apply the scaling ansatz, one must know the characteristic length at time r. Often R(t) is estimated 
from the structure factor, however we have found an analytic expression that works well. Indeed, one advantage of our 
numerical approach, as compared to Monte Carlo or cell dynamical simulations, is that we can test theories without 
using any free parameters. In particular, Ohta, Jasnow, and Kawasaki Jig] found that the time dependence of the 
domain size obeyed R(t) = ^/ApaT, with p c i = (d — l)/d. Furthermore, they found that the structure factor scaled, 
and they gave an explicit form for the scaling function Fi(t). In comparison, the theory of Kawasaki, Yalabik, and 
Gunton gives the same form of but the factor pd does not appear in their result for R(t). Our present 

simulations agree with the amplitude given by Ohta et ai, and so we choose a scaled time t given by 

t(k,T) - [kR{T)f = A Pd k 2 T . (2.8) 



III. TWO-TIME CORRELATION FUNCTIONS 



A quantity to which experiments give ready access is the fluctuation in the speckle intensity as a function of time. 
The relationship between an individual speckle at two different times T\ and T2 is, on average, described by the 
intensity covariance, 
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C0V k (k, 71,7a) = (/(k,7l)I(k,T 2 )) - (7(k,T 1 )^/(k ! T 2 )) . (3.1) 

For random systems, the covariance is maximum in the equal-time limit, t\ = 72, and as the two measurement times 
become widely separated, the values of the intensity become stochastically independent and the covariance decays to 
zero. For this relaxational system, negative values are not expected. The scaling ansatz extended to this situation 
allows collapse of the covariance at different (k, 71, T2) by 

k 2d Cov u (k, 7-1,730 = Cov(ti,t 2 ) (3.2) 

where Cov(ii,i2) is the scaling function for the covariance. 

The simulated scattering intensities were analyzed in the following way. In Eq. ([T^), the ensemble average over 
initial conditions is scaled to the universal form. To make estimates of Cov(ti,£ 2 ) from the simulations, we changed 
the order of the averaging and used the single-simulation averages Mi(t) and M 2 (ti, t 2 ). M\(t) is the scaled intensity, 
k d I(k, t), averaged over all pairs of (k, r) that map onto t, and M 2 (ti,t 2 ) is the scaled product of the intensities at 
two different times, k 2d I(k, ri)/(k, 72), averaged over all triples (k, ri,r 2 ) that map onto (iijta). Due to the large 
amount of data involved in the present study, samples for Afi(ii), Mi (£2) and M 2 (ti,t 2 ) were accumulated in a two- 
dimensional structure of bins organized by a pair of variables related to (ti , £2) as described below. After accumulation 
into the bins (ii,f 2 ), averaging over independent runs, i.e., over different initial conditions, was performed to further 
improve our statistics. These averages were then used to find the scaling function Cov(ii, t 2 ). It should be noted that 
all the quantities used here are also readily obtained experimentally, except for two unknown proportionality constants 
that depend on details of the experimental system. One occurs in the equation for the scaled time, corresponding 
to Apd in Eq. (^^). The other is the proportionality constant between the measured scattering intensity and the 



squared- norm of the Fourier transform of the order parameter, which we have ignored in Eq. (1.2). Neither of these 
should be a barrier to comparing our results to experiments. 

The normalized analog of the covariance is the correlation function |2C 

(l(k,7l)/(k,T 2 )) - (j(k,7l))(/(k,7- 2 V 

Corr(k,7i,7- 2 ) - N 2 / 2 ■ (3 - 3) 



/ 2 (k,ri)) - <;i(k,7i)) yj {P(k,r 2 )) - (/(k,r 2 ) 

Corr(k, t, t) is unity by construction, which removes the equal-time variations in the covariance as r changes. Using 
the definitions of Mi and M 2 above, the scaled version can be expressed as 



/ M 2 (ii,i 2 ))-(Mi(i 1 )\(M 1 (t 2 )\ 
Corr(ti,t 2 )= : v 1 \ Q 1 _ , (3.4) 



^M 2 (t u h) - <Afi(ii)) 2 ^M 2 (t 2) t 2 ) - (Mi(t 2 )) 2 



where (• ■ •) denotes averaging over initial conditions as before. The contour plot of Corr in the (ti, f 2 )-plane, Fig. 4, 
shows how the correlations in individual speckle intensities decay. In this figure, the line Corr(£i,i 2 ) = 1 extends 
along the diagonal. Moving away from that line, contours at values of 0.7, 0.2, 0.07, and 0.02 are shown. The scatter 
in the data is apparent for the last contour and becomes dominant for values of the correlation less than that. The 
striking feature of this figure is that the correlations increase in a nontrivial way as the phase-ordering continues. 
Thus, normalization by the time-dependent intensity is not sufficient to convert the speckle intensity to a stationary 
time series. 

A more natural set of variables for studying this effect is t = (ti + t 2 )/2 and 6t = \t% — ti\. A constant value of i 
corresponds to a line perpendicular to the t\ = t 2 diagonal, while 5t measures the distance (in units of scaled time) 
away from the diagonal. The symmetry under exchange of t\ and t 2 is retained. These variables were used for the 
binning of simulation results that produced Fig. 4, with data grouped into 250 equally wide series for < t < 1000; 
each t series having 100 equal-width bins with < St < 2i. 

The characteristic time difference, St c , required for the scaled intensity covariance, Cov(ii,£ 2 ), to decay to half its 
maximum value can be found as a function of i. (The normalized correlation function, Corr(ii,£ 2 ), does not decay 
to 1/2 for small values of t .) Here, results for Mi and M2 for < i < 2 (in 10 series, with 25 bins within each 
series) were collected in addition to those previously mentioned. The value of St c for each i was determined by linear 
interpolation, and the dependence of 8t c on i is presented in log-log form in Fig. 5. In this figure, two asymptotic 
limits giving different algebraic relationships are obvious. For small values of i the relationship is linear, with the 
exponent determined by least squares being 1.00 ± 0.01 for i < 1. At large values of i, a least-squares fit for t > 200 
gives an exponent of 0.49 ± 0.02. The estimates of error here are from obtaining the exponents from different ranges 
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of i in the simulation data; statistical error is an order of magnitude smaller than these estimates. The exponents 
we obtain are in good accord with theories discussed in Sec. below, which give them to be 1 and 1/2, respectively. 
However, the connection to theory requires some further justification, given in the following section. 



IV. CORRELATIONS IN THE SCATTERING MATERIAL 



While experiments can measure time correlations readily through the covariance or the correlation of intensities, Cov 
or Corr, theories to date have made use of the two-point, two-time order-parameter correlation function, C(r, n, r 2 ) = 
(?/>(0, n)ip(r, T2)). In this section we argue for, and empirically demonstrate, equality between the intensity covariance 
for k 7^ 0, which involves fourth moments of the order parameter, and the square of the spatial Fourier transform of 
the two-time order-parameter correlation function in our model. 

The relationship is exact when t/>(k) is a joint Gaussian random number, with its real and imaginary parts in- 
dependent and Gaussian, as we will show below. However, establishing ip(k) to be approximately Gaussian in the 
present case is not trivial. Gaussian variables are a natural consequence of the central-limit theorem, which requires a 
large number of uncorrelated contributions to the variable (with some restrictions on the properties of the individual 
contributions). For example, in a disordered system in equilibrium, correlations exist only on the scale of a small 
length £, so a system of edge-length L consists of on the order of (L/£) d independent, uncorrelated parts. Then the 
central-limit theorem applies, and ?/>(k) is a complex, Gaussian variable. For an ordered system in equilibrium, this 
argument applies to the fluctuations around the ordered state. The present nonequilibrium situation has some rough 
analogies to a disordered equilibrium state. At a given time, the average domain size is R(t), and so the number of 
independent parts at a given time is approximately (i/i?(r)) d , which can be large. It therefore seems reasonable to 
expect the central-limit theorem to apply, and indeed, we find empirically that ip()i) is Gaussian. However, unlike a 
disordered system, the distribution of domains of different sizes is broad in this case due to the initial long-wavelength 
instability and, furthermore, it is clear that domains interact as they grow. The degree to which this correlation 
is important is a nontrivial issue; below, we test it numerically and find these correlations to be negligible for the 
two-point quantities of interest in this work. 

If ip(k) is Gaussian, it is straightforward to relate the intensity covariance Covt to the order-parameter correlation 
function C. Wick's theorem can be used to decompose the intensity-intensity average as 

(7(k, Tl )7(k,r 2 )) = (v3(k,T 1 )^(k,r 1 )^(k,r 2 )^(k,r 2 )) (4.1) 

(^(k,r 1 )^(k,r 1 ))(^(k,r 2 )V-*(k,r 2 )) (4.2) 

+ (^(k,r 1 )^(k,r 2 ))(^(k,T 1 )V3(k,r 2 )) 

+ (Vi(k, n)^(k, T 2 ))(^*(k, Ti)^*(k, r 2 )) 
= (l + (5k,o)5 2 (k,r 1 ,r 2 ) + 5(k,T 1 )5(k,r 2 ) . (4.3) 

Here S*(k, t\, t 2 ) is the two-time structure factor corresponding to the two-point, two-time order-parameter correlation 
function 



5(k,n,r 2 ) = ^(k,Ti)V(k,7^ (4.4) 
= J dre^ r C(r,T U T 2 ) . (4.5) 



For k 7^ 0, Eq. (4.3) can be rewritten as 

C0V k (k,T 1 ,T 2 ) = 5 2 (k,T 1 ,T 2 ), (4.6) 

which equates the speckle intensity covariance with the square of the two-time structure factor of the system. Finally, 
the scaling ansatz defines a universal two-time form 

F 2 (t 1 ,t 2 ) = k d S(k,T 1 ,T2) (4.7) 

and 

Cov(t 1) t 2 ) = F 2 2 (t 1 ,t 2 ). (4.8) 
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Our numerical tests show that this equality holds well in these simulations for k ^ 0. In computer simulations, 



unlike scattering experiments, the two-time structure factor can be found directly using Eq. (i.4). Generally, the 
product ^(k,n)^*(k, t%) is a complex number, but the mean value of the imaginary part is zero, so the two-time 
structure factor is real valued. In our simulation data, the imaginary part of S(k,Ti,T2) is found to be zero, within 
our accuracy. Results for the real part at t\ — 25 and r 2 = 50 are presented on a log-log scale in Fig. 6. The direct 
measurements of the two-time structure factor and the square-root of the intensity covariance agree quite well, except 
at very large values of k, where lattice effects are important. 

Another consequence of our pro pose d decoupling, leading to the relationship between Covk and S, can be tested 



by simulation. When t\ — T2, Eq. (4.6) is simply 

(7 2 (k,r))-(/(k,r)> 2 = (/(k,r)) 2 , (4.9) 



and scaling can be applied to give F\ (t) = F% (t, t) — J Cov(t, t). These are compared in Fig. 7, w here the agreement 
is clearly quite good. As an aside, since the variance does not depend on the system size, Eq. (19) demonstrates 
that J(k, r) is a spatially non-self-averaging quantity pi] . Because of this, increasing the system size L will not 
improve the estimate of 5*(k, t) given by a specific number of speckles. However, this is compensated by the fact that 
more independent speckles are available in a given range of k for each trial. Numerical results obtained by Shinozaki 
and Oono in a study of spinodal decomposition in three dimensions using the cell-dynamical method appear to be 
consi sten t with this result |22fl . We al so n ote that the normalization of the two-time intensity correlation function, 



Eq. fl3.3|) , can be simplified using Eq. ( |4.9| ) 



Equation ( |4.9| ) is a property of exponentially distributed variables, and to the extent that ^>(k, r) is a Gaussian 
random number, /(k, r) will be an exponentially distributed random number. That is, the probability density for the 
normalized intensities, s(k, r) = /(k, r)/5(k, t), should satisfy 

P(s) =exp(-s) , (4.10) 

independent of (k, r). The probability density P(s) is normalized and has unit mean and standard deviation. Since 
P(s) is identical for all values of (k, r), only one density function needs to be constructed. The results for all 
0.024 < k(k) < 0.75 are presented for two times in log-linear form in Fig. 8. The histogram for s is constructed 
with a bin size of 0.1, and the normalized intensity is found using the circular average of /(k, r) from the same trial. 
The lower bound on k is chosen such that at least 20 speckles contribute to this circular average. This histogram is 
accumulated over all 100 trials and then normalized; for each histogram on the order of 10 6 samples are available. 
The solid line is the expected density, P(s) = exp(— s). For the earlier time the probability density is exponential for 
all s. At the latest simulation time the density is exponential only for s < 5, with a higher than expected occurrence 
of speckles brighter than this. Still, less than 0.5% of the points lie this far into the tail, and the deviation is inferred 
to be a finite-size effect because it becomes more pronounced as our simulation continues. This agrees with earlier 
simulations fTq ], in which increasing the system size eliminated the deviation. We believe this effect is an artifact of 
the periodic boundary conditions which allow stabilization of slab- like domain patterns pq ], producing "frozen- in" 
interference patterns that remain bright while the average intensity steadily decreases. No k dependence is found for 
P(s). 

V. TWO-TIME THEORIES 

Two theoretical predictions for the two-point, two-time order-parameter correlation function exist, which will be 
described and then compared to our simulation results. The first is an analytic theory, due to Yeung and Jasnow p4j , 
which is an extension of the analysis by Ohta, Jasnow and Kawasaki fjgj , p5[| . The second theory, developed by Liu 
and Mazenko ]26|], is numerical. Both theories produce approximations for the two-point, two-time order parameter 
correlation function, C(r, Ti,T2), in the asymptotic limit. These then give the structure factor through the general 
relation for the Fourier transform of a spherically symmetric function /(r) in d dimensions, 

J dre lkr f(r) = (2Tr)%k- d duui J&_ x {u) f{u/k) , (5.1) 

where J v is a Bessel function of the first kind of order v. 

In the scaling regime, where the domain-wall thickness can be ignored, the Yeung- Jasnow correlation function is 
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/ 2R( n )R(r 2 ) 



V^n) 2 + R( T2 y 



exp 



R(nf + R{T 2 f 



CYj(r,Tx,T 2 ) = -arcsin 

TT 

With the previously defined variables, St = \t\ — t 2 \ and i = (t± + t 2 )/2, this yields 

d 

2 ^ roc ^ / / St\^\ 4 

F 2t yj(St,i) = — (27r)' J du u 2 Jd„ 1 (u) arcsin I 1 — J j ex P 
Expansion of arcsin(x) about x = 0, followed by term-wise integration, gives 

d(2j+l) 



2i 



F a ,Yj(«t,f) = -(27r)*t«y; 



(2i)! 



=J 2W(7!)2(2j + l)- 



■ exp 



2(2j + 1) 



(5.2) 



(5.3) 



(5.4) 



This infinite series is convergent for all physical values of St and i; however, its terms are in general nonmonotonic 
in j. Only as one or both of the scaled times t\ and t 2 approaches zero, is the series well approximated by the j—0 
term. The asymptotic early-time form, 



F 2 MSt,t) ~ fa)! 1-(|) 2 *<exp(-|), 



(5.5a) 



is therefore a good approximation to Eq. (5.3) only in the very restricted region, 

2 



< 1- 



<gC 6 d 3 exp 



-2t 
~3d 



(5.5b) 



near the t\ and t 2 axes or the origin. 

A more useful, analytical result is obtained in the limit of large t and small St. Then, the largest t erm s in the series 
occur in a relatively wide range of j near i/2(d + 3). For large t, the exponential factors in Eq. (5.4) suppress the 
small- j terms. The series can then be converted to an integral, and the factorials can be approximated by Stirling's 
formula to give 



F 2 ,Yj(St,t) « 7T _5_ 4 2 t 2 



du- a- - 1 1 — [ 

2t 



2\ 8^; 



exp (— w) 



(5.6) 



The identity lim m ^ 00 (l — x/m) m = exp(— x) is used to obtain the explicitly integrable form, 

z 2 d 



d-3 d+1 - 1 
F 2 ,Yj(Z,t) W 7T 2 4 2 < 2 



du> u> 2 exp 



32u> 



(5.7) 



wher e z = St/vt. We note that the full Yeung-Jasnow result, Eq. (5.3), as well as the early-time approximation, 
Eq. ( 5.5a ), depends on St only through the scaling combination St/t. However, in the asymptotic late-time approx- 



imation, Eq. (5/7), the natural scaling combination is z — St/y i. A s we shall see below, this analytical result is in 
excellent agreement with our numerical simulations. Equation ( |5.7| ) is readily integrated to yield the explicit large-f 
asymptotic scaling function, 



t*F 2 , Y j(z,t) w TT— T 



2 z\ 



K d+i z\ 



(5.8) 



where K n is a modified Bessel function of the second kind. The right-hand side of this equation depends only on z. 
In this large-t limit, the asymptotic forms of this scaling function with respect to z are 



i*F 2tYJ (z,i) 



16 d-l' 
Z\ 



for z«l 

^2^(z v /f) f exp(-z^l) for 1 « z « 2y/t 



(5.9) 
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The first line of this equation is exact for z — 0, where it gives the asymptotic Porod-tail limit of the Ohta-Jasnow- 
Kawasaki result for t^Fify). Note that, in the limit St = 2t, the second line gives a nonzero value, in contrast to the 
proper resu lt given in Eq. (5.5a). When z is on the order of yt, i^.Yj diverges from the small-z approximation given 



inEq. (kg). 

A nonrigorous scaling argument ]2q] suggests that the order-parameter correlation function should depend on 
wavevector and time through fc|-R(ri) — R(t2)\ oc St/t^~ n \ when St <C t. The scaling variable z obtained above 
is consistent with this since n — 1/2 for Model A. In fact, the z — <5t/f( 1 ~™" ) result is obtained in the asymptotic 
large-i limit if one repeats the above calculation for general n, considering the Yeung-Jasnow form for the correlation 
function, Eq. (^2), simply as an integrable approximation valid for small r. 



The second theory for two-time correlations in Model A is due to Liu and Mazenko [[26|. It is an extension of a 
theory developed by Mazenko J27|] to predict the universal part of the two-point, one-time scaled order-parameter 
correlation function C(r/R(r)). The heart of the Liu-Mazenko theory is the scaling ansatz 

C(r,ri,r 2 ) - C LM (r/R(n), Wn) (5.10) 

and the partial differential equation 

dCh ^ ^ = V|Clm(x, /) + 2x ■ V a C L M (x, r') + ^ tan (|c LM (x, r')) . (5.11) 



In Eq. ( [5.11 ), Clm is considered a function of the rescaled variables x = t/{2^/t) and At' = 1h(t2/ti). Note that 



these definitions are in terms of the physical quantities, but the connections to the rescaled variables used in this 
work are straightforward. We have found numerical solutions in terms of Liu and Mazenko's units, which we then 
converted so that all results shown here are in terms of the rescaling given in Section ||. The Mazenko theory for 
one-time correlations pTJ serves to provide input into the Liu-Mazenko model through the scaled order-parameter 
correlation function Cm(x) = Clm(^, t'=0) and the numerically obtained eigenvalue /i*. Since the system is assumed 



isotropic, Eq. (5.11) can be reduced to a single partial differential equation in terms of a radial distance x and the 
logarithmic time ratio r'. 

The scaled order-parameter c orrel ation function in the Liu-Mazenko theory is normalized to Clm (0,0) = 1, which 



causes the tangent term in Eq. (5.11) to diverge. Use of the transformation G — sin (7tClm/2) leads to the small-time 
solution 

C LM (0y) = ^arcsin(exp (- ^Tf') ) ■ (5-12) 

This solution is only needed to avoid numerical difficulties at x = for the first time increment. Aside from this, 
Eq. ( p. 11 ) can be solved numerically using a finite differencing scheme that is implicit with respect to the derivatives, 



but evaluates the tangent term explicitly. Using Ax = 0.01/i* and At' = 10 we have reproduced Fig. 1 of 
Ref. p(|. In addition, the exponent for the autocorrelation Clm(0,t') (which will be discussed later) is recovered. 
The Liu-Mazenko results can be compared to the theory presented here by taking t to be a parameter and noting 
that St — 2t tanh (2t ! ) . The two-time structure factor predicted by the Liu-Mazenko theory, found using the Fourier 



transform in Eq. (5.1), is 



F 2}h u(St(i,T'),t) = (2tt)* [ —-JL—— ) x (5.13) 

\ p d (l + exp (-4r')) / 




dxx " J ^ 1 x '! p d (i + Z(-±T')) ) CMx > T ' ] ' 

where the scaling function Clm(^, t') itself depends on d. 

The two-point, two-time order-parameter correlation functions predicted by the two theories can be compared 
directly with our simulation, without adjustable parameters. Both theories describe the data well in some instances, 
and poorly in others. Our results are presented in Fig. 9 for t\ = 100 and ti = 200. Data from our simulation are 
circularly averaged with bins of width one in the rescaled distance units, then averaged over 80 trials. The agreement 
between the Yeung-Jasnow theory and our simulation is quite good for this choice of times, while the Liu-Mazenko 
result is only qualitatively correct. 
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For larger separations in time, the Yeung-Jasnow theory does not work as well, most noticeably in predicting the 
autocorrelation lunction, 

A(Ti,7a) = <V(r,Ti)V>(r,7*)) , (5.14) 

which is equivalent to C(r = 0,ti,T2). Fisher and Huse pj| have argued that for t 2 3> t\ the autocorrelation obeys 
the power-law 

A( n ,T 2 )~[R( n )/R(T 2 )] X (5.15) 

for late times. (Note that, since R oc t™, the variables used in this paper give In [R(ti) / R(t 2 )] = — 2rttanh~ 1 St/2i.) 
Fisher and Huse give physical arguments for d/2 < A < d. They point out that A = 1 for the d = 1 Glauber model 
and conjecture that A = 5/4 for the two-dimensional spin- flip Ising model. In the limit t 2 3> ti, the Yeung-Jasnow 



theory gives A = d/2, as is easily seen by setting r — in Eq. (5.2). The Liu-Mazenko theory yields A « 1.2887 and 
1.6726 for two and three dimensions, respectively p6| . Cell-dynamical simulations performed by Liu and Mazenko 
p6[ gave A = 1.246 ± 0.02 for d — 2. A recent experiment |3(J found A = 1.246 ± 0.079 for a two-dimensional nematic 
liquid crystal using video techniques. For our simulations we measured A(ti,t 2 ) for several values of n. The results 
are presented on a log- log scale versus t 2 /t\ in Fig. 10. [In this measurement we did not employ the nonlinear 
transformation tp — ► sign(V>). This accounts for the fact that A{tx,Tx) < 1, but does not otherwise seem to affect 
the results.] In the figure, the data appear to support a power-law decay at the latest time, and the exponent found 
by fitting the 48 points t 2 > 1800 for t\ = 20 is A ps 1.24, which is in good agreement with the experiment and the 
simulations by Liu and Mazenko. However, the local effective value of A, A e ff = 2<ilnyl/<iln(Ti/T2), which is shown 
vs. t 2 /t\ in the inset in Fig. 10, does not show a clear convergence to an asymptotic limit, especially for larger values 
of t\ . Here A c g is obtained as a three-point finite-difference estimate around t 2 , but estimates that smooth the data 
over wider intervals yield similarly irregular results. The estimates of A obtained from our present simulations are 
thus somewhat uncertain. However, the Yeung-Jasnow prediction, A = 1, is clearly violated. 

The predictions of the two theories for the two-time structure factor are compared to our simulation data in Fig. 6. 
The Liu-Mazenko theory is in semi-quantitative agreement with our simulations at small k, but falls off much more 
rapidly at large k. In addition the shoulder present both in our simulation result and in the Yeung-Jasnow theory 
is absent in the Liu-Mazenko result. The Yeung-Jasnow theory agrees much better with our simulation, although it 
does not fall off as fast at large k, and its overestimation around the shoulder is seen consistently throughout our 
simulations. 



The prediction for the characteristic scaled time separation, St c , defined in Sec. [II, for both theories is compared 
to our simulation results in Fig. 5. All three agree that St c = Ddi for small t. Least-squares fits give D 2 = 1.36 ± 0.02 
for our simulation, D 2 = 1.325 ± 0.002 and D 3 = 1.633 ± 0.006 for Yeung-Jasnow, and D 2 = 1.00 ± 0.02 for Liu- 
Mazenko. At large i, the St c = Ddi 1 / 2 behavior of the Yeung-Jasnow approach results naturally from Eq. fl5.8| ). This 
equation can be solved numerically to find D 2 ps 2.156 and D 3 2.187 in two and three dimensions, respectively Our 
simulation results give D 2 = 2.12±0.01. The value off separating the two scaling behaviors corresponds quantitatively 
to the value for which the first term no longer dominates the series expansion of the Yeung-Jasnow result, given by 
the relation ( |5.5b| ). The agreement between the Yeung-Jasnow theory and our simulation is quite remarkable. On the 
other hand, for t not small the Liu-Mazenko theory crosses over to a region where 5t c is independent of t and then 
into another linear region; neither of these relationships is seen in our simulations. 

The analytic expression for the uni versal form of the two-time structure factor deduced from the Yeung-Jasnow 
theory for large i is given in Eq. ( |5.8| ). It is tested for several values of t in Fig. 11, which shows good collapse of 
the simulation data for tCov(<5t, i) in terms of the scaling variable z — St/t 1 / 2 . The data also agree quite well with 
the Yeung-Jasnow scaling function, iF 2 yj f° r z < 5. Indeed, the slow quadratic decay of correlations near z = is 
another signature of persistence in the phase-ordering system. In contrast, Brownian fluctuations give exponential 
decay from St = 0. The agreement with the Yeung-Jasnow theory is remarkable since the scaling of the simulation 
data uses an analytic expression for the characteristic length, and no adjustable parameters are employed. 



VI. CONCLUSIONS 



Using both numerical and analytic methods, we have investigated time-time correlations in the scattering intensity 
for a two-dimensional system undergoing an order-disorder transition. The correlations are found to obey scaling in 
terms of the variables 5t = \t 2 ~ ti\ and i = (t\ + t 2 )/2. In the large i limit, the correlation data collapse onto a 
universal curve which is a function only of 5t/Vt. 

We argue for, and establish numerically, an exponential distribution for the scattering intensity, Eq. ( 4.10| ), and 
equality between the scattering intensity covariance and the square of the two-time structure factor of the order 
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parameter, Eq. (4.6), for k ^ 0. We use this equality to test theories for the two-time structure factor due to Yeung 
and Jasnow, and Liu and Mazenko. Both theories describe the data well in some instances, and poorly in other 
cases. The Yeung- Jasnow theory is very similar to our simulation results, so long as St is not too large. For T2 Ti, 
the Liu-Mazenko theory gives a better estimate for the autocorrelation exponent A. For large t, however, the Liu- 
Mazenko theory does not show the same scaling as our simulation results, where the Yeung- Jasnow theory compares 
quantitatively well. 

Our numerical simulations indicate that a definitive experimental treatment of time-correlations during an order- 
disorder transition is possible by intensity-correlation spectrometry of scattering speckle. Analysis of experimental 
correlation data should be similar to the procedures discussed for the simulation data in Sec. III. For non-conserved 



systems, the experimental scaling function should be well approximated by Eq. (5.5a), for small fc, with one adjustable 
para met er for each axis. With a similar adjustable parameter scheme, the scaling function should be described by 
Eq. (5.8) for data in the Porod tail. 



Finally, we expect that the equality between the intensity covariance and the squared two-time structure factor 
also occurs in other phase-ordering systems; in particular, we expect that it occurs for conserved systems. That 
would allow the experimental study of, for example, time correlations in binary alloys undergoing phase separation 
by spinodal decomposition, which are representative of Model B. Indeed, preliminary numerical work we have done 
indicates this. Experiments on such systems would be of considerable value. 
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FIG. 1. A typical configuration of domains, taken from one of the simulations reported here. Here all systems are 1024 x 1024, 
and this picture is for the latest simulation time, r = 2000. 




FIG. 2. Example of the speckled scattering intensity; this pattern corresponds to the domain structure shown in Fig. 1. The 
intensity is shown on a logarithmic scale with darker shades indicating brighter speckles. This is a 200 x 200 section from 
the full 1024 x 1024 pattern with the k = origin at the center of the figure. Speckles do not shift in k-space, but their 
intensities fluctuate strongly around the (k, r)-dependent average value. The rays are present in individual patterns, but are 
not correlated between trials. 
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FIG. 3. Time evolution of the scattering intensity at one wavevector for one quench to zero temperature. The intensity 
has been normalized by the time-dependent structure factor determined over all 100 simulations considered here. The dotted 
line is a synthetic "Brownian" function constructed to have an exponential single-time probability density and an exponential 
two-time covariance with a characteristic time corresponding to that of the simulated intensity. Compared to the "Brownian" 
function, the persistence of the scattering intensity is apparent [JL 3|Jl4[] . 
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FIG. 4. Contour plot of the scaled two-time intensity correlation function, Corr(ti,i2)- The correlation at t\ = <2 is unity by 
construction. The contours moving away from this diagonal are at 0.7, 0.2, 0.07, and 0.02. The figure shows that the speckle 
intensity stays correlated for larger values of St = |i 2 — ti| as i — (ti + t-i)j2 increases. 
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FIG. 5. The characteristic decay time difference, 5t c , as a function of i = (ti +t2)/2. Power-law behavior is seen at small and 
large values of t. Least-squares fits to the simulation data yield a linear relationship at small i and St c ~ i a with a « 1/2 at large 
i. The heavy solid line is the analytic prediction of the Yeung-Jasnow theory. The broken line is its small-f limit St c ~ 1.325i, 
and the li ght so lid line is the large-t approximation St c w 2.156%/T. The the t associated with the change in behavior agrees well 
with Eq. (5.5b). The results of the Liu-Mazenko theory are represented as circles connected by a dot-dashed line; least-squares 



fits show them to be nearly linear at both small and large t. 




FIG. 6. The two-time structure factor for ti = 25 and T2 = 50. The simulation results are found directly as S(k,Ti,T2) 
(squares), and indirectly as v/ Covk(fc, ti, t-i) (diamonds). The Liu-Mazenko theory agrees semi-quantitatively at small k. The 
Yeung-Jasnow theory agrees much better, even reproducing the second shoulder qualitatively. 
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FIG. 7. Estimates of the scaling function for the structure factor, Fi(t), obtained from the simulation data using the average 
scaled scattering intensity, (Mi(t)) (circles), and the square-root of its variance, y / Cov(i,t) (diamonds). The agreement between 
the two measurements is a direct consequence of the intensity being an exponentially distributed random variable. The forms 
obtained from the Yeung-Jasnow (solid line) and Mazenko (dashed line) theories are also included. 
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FIG. 9. The two-time correlation function C(r, T\,T2) with n = 100 and T2 = 200 presented on a semi-log scale. The 
characteristic length R(t 2 ) is determined analytically as described in the text. The agreement between the simulation (circles) 
and the Yeung-Jasnow theory (solid line) is quite good for this choice of parameters; the Liu-Mazenko theory (dashed line) is 
noticeably different. The agreement between simulations and Yeung-Jasnow is not as good for larger time separations. 




FIG. 10. The autocorrelation function A(t\,T2)- A least-squares fit for T2 > 1800 for the n = 20 data gives the exponent 
A w 1.24. Fits for later ri give progressively smaller values, but are still larger than 1. The latest time value is approached 
from below, and the asymptotic value may not be obtained for the simulation times considered here. For two dimensions, the 
Liu-Mazenko value of A is found numerically to be approximately 1.2887, while a recent experiment gives 1.246 ± 0.079. The 
Yeung-Jasnow prediction is unity. The local effective exponent, A e ff, estimated from a three-point finite difference appear in 
the inset. They characterize the uncertainty in our estimate of A, but show that the Yeung-Jasnow prediction is violated. 
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FIG. 11. Plot of the intensity-covariance scaling function, tCov(St,t), versus z = 5t/vt for different values of t. The data 
collapse onto a single curve for a large range of i. T he solid curve is the corresponding analytic scaling function, tF^YjiSt^i), 



predicted from the Yeung-Jasnow result in Eq. (5.8). It agrees quite well with the simulation data. 
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